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O • Abstract. We consider 2-state quantum walks (QWs) on the line, which are defined by two ma- 

^__> ■ 

trices. One of the matrices operates the walk at only half-time. In the usual QWs starting from the 
origin, localization does not occur at all. However, our walk can be localized around the origin. In 
this paper, we present two limit theorems, that is, one is a stationary distribution and the other is a 
^! ■ convergence theorem in distribution. 
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1. INTRODUCTION 



^ ■ The quantum walks (QWs) are considered as the quantum counterparts of a random 



walk. The discrete-time QWs introduced by Aharonov et al. [1] and Meyer [2] were 
intensively studied in Ambainis et al. [3]. On the other hand, the continuous-time models 
were defined by Farhi and Gutmann [4]. We concentrate on the discrete-time 2-state 
QWs in this work. In relation with quantum computer, QWs are often investigated. Since 
the quantum walker has ballistic behavior in probability distribution, quantum searches 
which are more fruitful than the classical ones can be designed by QWs [5, 6, 7, 8]. 



In other research field, applications of the walks were also discussed [9, 10, 11]. The 
asymptotic behaviors of the walks have been analyzed, and limit theorems have been 
obtained. For example, Konno [12, 13] calculated the limit distribution of the usual 2- 
state walk (i.e., time-homogeneous model) on the line Z = { 0, ± 1 , ±2, . . . } which had a 
simple density function. In the present paper, we consider a time-inhomogeneous 2-state 
quantum walk (QW) determined by two unitary matrices on the line. There are some 
results for time-inhomogeneous models [14, 15, 16]. Banuls et al. [14] and Romanelli 
[16] insisted localization for their models, respectively. Machida and Konno [15] gave 
convergence theorems to a 2-period QW with two unitary matrices and two special 
time-dependent walks. We focus on calculation of the limit distribution which can be 
expressed by a combination of density functions. QWs whose limit distributions had 
similar structures to ours were treated in [17, 18, 19]. Brun et al. [17] introduced and 
discussed the multi-state walks with the tensor-product representations. A limit theorem 
for the multi-state walk with a special initial state was obtained in Segawa and Konno 
[19]. Miyazaki et al. [18] also computed the limit distribution of the walk defined by 
the Wigner formula of rotation matrices. If the number of states is more than four in 
these multi-state models, the limit distributions are depicted as a combination of density 
functions. 

This paper is organized as follows. In Sect. 2, we define our walk. Sect. 3 is devoted 
to calculation for the Fourier transform of amplitude. We present limit distributions as 
our main result in Sect. 4. A detail of the computation for the limit distributions is given 
in Appendix. By using the Fourier analysis introduced and studied by Grimmett et al. 
[20], we obtain the limit distribution. Summary is given in the final section. 



2. DEFINITION OF OUR 2-STATE QWS 

In this section we define our 2-state QWs on the line. Let |jc) (x G Z) be an infinite 
components vector which denotes the position of the walker. Here, x-ih. component of 
\x) is 1 and the other is 0. Let | £ C 2 be the amplitude of the walker at position 

x at time t e {0, 1, 2, . . .}, where C is the set of complex numbers. The walk at time t is 
expressed by 

w = Ew®i^w>. (i) 

The time evolution of our walk is constructed by the following two unitary matrices: 



u = 


cos 6 


sin0 




c s 


,H = 


1 




sin0 


— COS0 




s —c 




-1 



where c = cos , s = sin and € [0, 2tt) (0 ^ 0, f , tt, ^) . Moreover, we introduce four 
matrices: 





c s 









1 







p = 





,2 = 


5 — C 


,P1 = 





,<2i = 


-1 



Then the evolution is determined as follows: 
For? = t+1,2(t+1),...,wi(t+1), 



V^(x))=Pi|^_i(x+1)) + (2i]^_i(x-1)), (4) 



and for T+ 1,2(t+ 1), . . . ,m(z+ 1), 



|^(x))=P|^-i(x+l)> + Q|^-i(^-l)), (5) 

with m, T e {1,2, . . .}. Note that P + Q = U and Pi + Q\ = H. The probability that the 
quantum walker X t is at position x at time t, ¥(X t = x), is defined by 

P(X f =x) = (^ f (x)|^(x)). (6) 

The Fourier transform \^ t (k)) (k G [-7C, it)) of | %{x)) is given by 



l*r(*)> = E«~ to lv*(*)>- (7) 



By the inverse Fourier transform, we have 



/% Air 
/ kx \*t(k))-. (8) 



From Eqs. (4), (5) and (7), the time evolution of |*P f (fc)) becomes 



(? = T+l,2(r+l),...,m(T+l)), 
= < (9) 

U(k)\%^(k)) (^t+1,2(t+1),...,hi(t+1)), 



where U(k) = R(k)U , H{k) = R(k)H and R(k) 



e ik 



In order to analyze the 



e~ ik 

walk after time m(r+ 1), we focus on the probability distribution especially at time 
(m + n)T + m (n e {1 , 2, . . .}). Then the Fourier transform of amplitude can be written as 
follows: 

\V {m+n) r +m (k)} = U(kr(H(k)U(k)T \%(k)) . (10) 

In Sect. 3, Eq. (10) will be expressed by the eigenvalues and eigenvectors of U(k). In 
the present paper, we take the initial state as 

T [a,[5} (x=0), 
IW)(-v)) - { (ID 
r [0,0] (x^O), 

where |a| 2 + |j8| 2 = 1 and T is the transposed operator. We should note that \^o(k)} = 
IVo(0)>. 

3. THE FOURIER TRANSFORM OF AMPLITUDE EXPRESSED 
BY THE EIGENVALUES AND EIGENVECTORS OF U(k) 

In this section we rewrite | v I / ( m+ „) 1 : +m (x)) in Eq. (10) by using the eigenvalues and 
eigenvectors of U (k) whose formulae will be given in Appendix. At first, the Fourier 
transform at time m(r+ 1) can be described as follows: 



2 2 

I*»(t+i)(*)> = t (vj(k)\V m{T+l) (k)) \vj(k)) = £ a hm (k) \ Vj (k)) , (12) 

7=1 7=1 



where cij_ m (k) = (vj(k) 1^(^+1) (^)) an d \vj(k)) (j G {1,2}) are the normalized eigen- 
vectors corresponding to the eigenvalues Xj{k) of the unitary matrix U(k). In par- 
ticular, since the initial state |VbW) i s given by Eq. (11) in our model, we have 
ajfl(k) = (vj(k)\%(k)) = <Vj(ife)|Vb(0)>. From Eq. (12), we see that 



\v m{z+l) (k))=H(k)u(ky\v 

(m-l)(T+l) (*)> 

= ^a jhm - 1 (k)X h (kyH(k)\v h (k)) 

71=1 
2 2 

= E E ^h,m-i{k)x h {kyA nh {k)\v h {k)), 

71 = 1 72=1 



(13) 



where A hh (k) = {v h (k)\H{k)\v h (k)) (j h j 2 G {1,2}). Therefore we get 



ai, m (*) 
ai,m{k) 



h{kYA n {k) X 2 (kyA l2 (k) 

X l {kYA 2l {k) HkyA 22 (k) 

M(kyA n (k) X 2 (kYA 12 (k) 

x x {kyA 2l {k) x 2 (kyA 22 (k) 



a\, m -\{k) 
ai,m-\{k) 

aifl(k) 
a 2:0 (k) 



(14) 



Moreover, the coefficients a,j, m {k) of the eigenvectors \vj(k)) can be computed as fol- 
lows: 



ai,m(k) = 



z 2 (k)-zi(k) 



{{-lY{z2{k) m - l -zi{k) m - l )a hQ {k) 



a2,m(k) 



1 



zi{k) -z\{k) 



+ Ai (k) 'An(k) (Z2(k) m - zi (k) m )a lfi (k) 
+h(kyA 12 (k)(z 2 (k) m -zi(k) m )a 2fi (k)} , (15) 
{X x {kyA n {k){z2{k) m -zi{k) m )a hQ {k) 

+ (z2(kr +1 -z l (kr +l )a 2 ,o(k) 
-h{kyA n {k){z2{k) m -zi{k) m )a 2 ,o{k)} , (16) 



where Zj{k) (j G {1,2}) are the eigenvalues of the following matrix: 



XxikYAnik) X 2 (kyA l2 (k) 

h(kyA 21 (k) x 2 (kyA 22 (k) 



(17) 



that is, 



Zj (k) = v(k)-(-iyJv(ky+(-i)\ 



(18) 



and 



v(*) 



hikyAnW-i-iy^ikyAnik) 



(19) 



Hence, for a m j(k) given by Eqs. (15) and (16), Eq. (10) becomes 



\^( m+ n)r +m (k)) = (U(kyy |* m(T+1) (*)> = £ a hm (k)Xj(k) n * \ Vj (k)) . (20) 

7=1 



4. LIMIT DISTRIBUTIONS FOR THE WALK 



In this section the limit distributions for some cases, one of which will be computed 
in detail in Appendix, are presented. Firstly, we show the limit distribution for the usual 
walks which was given by Konno [12, 13]. Since the matrix H does not operate our walk 
till time t, we can obtain the convergence in distribution of the walk X x jx as the limit 
theorem for the usual 2-state walks. That is, for the walk, we have 



lim P ( — < x 

X—¥°° \ % 



fa{y\c) 



a 



l/3| 2 + 



dy, (21) 



where 



fn(x;c) 



7t(\ —X 2 )Vc 2 —X 2 



7 H4kl)M> 



(22) 



and Ia(x) = 1 if x e A, Ia(x) = if jc ^ A. Fig.l depicts the comparison between the 
limit density functions and the probability distributions at time 500 for usual walks with 
e = k/4. 



-1 -0.5 



-1 -0.5 0.5 1 



x t /t 



Xtft 



(a) |Yb(0)) = 1 [l/V2,i/V2] (b) |vb(0)> = J [1,0] 

FIGURE 1. Comparison between the limit density functions (thick line) and the probability distribu- 
tions at time t = 500 as T = 500 (thin line) with 9 = n/A. 



m = 1 case 



From now, we propose our results for the walks. At first, we concentrate on m = 1 
case and consider the probability distribution at time t = (n+ 1)t+ 1. The limit density 
functions as x — > °° are calculated for n = 1 case in which the limit distribution was 
computed by Machida [21] and n = 2, 3, . . . case, respectively. 
n — 1 case 

For n=l, Machida [21] obtained the limit distribution expressed by both a 5-function 
and a density function as follows: 



hmF (^I±l< x 



AMy)+f K (y;c)M l (y)dy 1 



(23) 



where 5q(x) is Dirac's 5-function at the origin and 



Ai 



i + kl' 



M l (x) = 



|a| 2 -|j8| 2 + 



l-^ z 



(24) 
(25) 



We should note that Ai doesn't depend on the initial state | yb(0)) = T [a, j8] . Assuming 
the initial state as | Vb(0)) = r [^7j , an d | Vb(0)) = r [ 1 , 0] , we get the limit densities 
and the probability distributions as shown by Fig. 2. 



-1 -0.5 

x t /t 



-1 -0.5 0.5 1 



x t /t 



(a) |y/b(0)> = r [l/v / 2,//v / 2] (b) |yb(0)> = r [l,0] 

FIGURE 2. Comparison between the limit density functions (thick line) and the probability distribu- 
tions at time t = 401 as z = 200 (thin line) with = n/A,m = !,« = !. 



n = 2,3,. . . case 

For n = 2, 3, . . ., the limit distribution is denoted by two density functions as follows: 



limP(5^±i<* 



/js- (y; c)Mi (y) + f K ( ^jy; c)B x ( ^J-y; «) rfy, (26) 
3 n— 1 n— 1 



where 



Bi(x',n) 



n+1 
n-l 



l + { |a| 2 -|/3| 2 + 



[afi + afi)s 



? 2 



(27) 



We show the symmetric and asymmetric distributions for n = 2 in Fig. 3. 



-1 -OS 05 

x t /t 



_ t». L — A J , 

1 -0.5 0.5 1 

Xt/t 



(a) \y (0))= T [l/V2,i/V2] (b) |yb(0)) = r [l,0] 



FIGURE 3. Comparison between the limit density functions (thick line) and the probability distribu- 
tions at time t = 601 as z = 200 (thin line) with = n/A,m = l,n = 2. 



m = 2 case 



Next we consider m = 2 case and focus on the probability distribution at time t = 
(n + 2)t + 2. In this case we compute the limit density functions as x — > °° for n = 1, 
n = 2 and n = 3,4, . . . cases, respectively. For n = 3,4, . . ., the calculation will be written 
in Appendix. 
n = 1 case 

For n = 1, the limit distribution has a combination of two density functions as follows: 



limPi 



< X 



/X 
My;c)M 2 ( y )+M3y;c)M 3 (3y)d y , (28) 
-oo 



where 



M 2 (x) 



l-< |a| z -|/3| 2 + 



2 (a/3 + a/3 > 



2 \ 2 

1 S 2 
1 ^X 



3 S 2 

M 3 (x) =— x 2 



2(a/3 + a/3>(l+s 2 ) s 2 2 



-3{ \a\ z -\p\ z + 



2 (a/3 + a/3)s I s 2 3 4(a/3 + a/3)s 3 x 



l-x 2 



(29) 



(30) 



When we take the initial state as | Vo(0)} = , -j=\ and | Vo(0)} = [1,0], the distri- 
butions are shown in Fig. 4 
n = 2 case 



For n = 2, we can describe the limit distribution with a 5-function and two density 



X t ft Xt/t 



(a) |^o(0)) = 7 [1/V2,i/V2] (b) |vb(0)) = '[1,0] 

FIGURE 4. Comparison between the limit density functions (thick line) and the probability distribu- 
tions at time t = 602 as z = 200 (thin line) with 6 = K/4,m = 2,n = 1. 



functions as follows: 



limP^^^<x)= / A 2 8o(y)+fK(y,c)M 2 (y)+f K (2y;c)B 2 (2y;2)dy, (31) 



where 



A2= c_M + (s*-c 2 )\ S \(\-\s\) 2 



n + 2s 2 2 
B2\x;n) = 



n cr 



2c 4 

l + i|a| 2 -|/3 



2 (a/3 + a/3>(3s 2 +r 



2 (afi + afi)s I s 2 3 4(a/3 + a/3)s 3 jc 



-4<^ |c^-|/3| 2 + 



1-x 2 



(32) 



(33) 



We should remark that A 2 is independent from the initial state | Vo(0)) = T [ a ? Fig- 5 
presents the limit density functions and the probability distributions. 



1 i. I a J 



-ML* 



-1 -0.5 



Xt/t 



-0.5 



(a) |vb(0)) = 7 [1/V2,//V2] (b) |vb(0)> = 7 [1,0] 

FIGURE 5. Comparison between the limit density functions (thick line) and the probability distribu- 
tions at time t = 802 as T = 200 (thin line) with 6 = %/A,m = 2,n = 2. 



n = 3,4,... case 

For n = 3,4, . . ., the limit distribution is expressed by three density functions as follows: 



limP(^±^<* 



/•■* n + 2 n + 2 n + 2 n + 2 

: / My;c)M 2 {y)+f K ( y;c)fl 2 ( >>;«)+/#( — ^j;c)5 3 ( — -y;n)dy, 

J-oo n n n — 2 n — 2 



where 



(34) 



B 3 (x;n) 



n + 2 
n — 2 



i i I? ioi2 (a/3 + a/3> 
1 + { \a\ - \p\ 2 + ^ J - — 



x 1 1 - 



(35) 



We draw the comparison between the limit density functions and the probability distri- 
butions for n = 3 in Fig. 6. 
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x t /t 



Xt/t 



(a) \ ¥o (0)) = T [l/V2,i/V2] (b) |y (0)> = r [l,0] 

FIGURE 6. Comparison between the limit density functions (thick line) and the probability distribu- 
tions at time t = 1002 as T = 200 (thin line) with 9 = n/4,m = 2,n = 3. 

We should note that the function fx ( ^^x; c)2?2 ( "jp*! ft ) can be finite at either x = -^-^ \ c \ 
or x = —jn^\c\. For example, if 9 = f and | Vb (0) ) = T [ i i0]' we see tnat 



lim f K ( H J '~^' x;c)B2( n x\n) = 0. 



(36) 



Fig. 7 shows typical examples of the distributions in the case when Eq. (36) is realized. 



-1 -0.5 \ 0.5 



X,/t 



UbLiL 



-1 -0.5 1 0.5 



x t /t 



(a) time t = 802 



(b) time? = 1002 



(T = 200,m = 2,n = 2) 



(T = 200,m = 2,n = 3) 



FIGURE 7. Comparison between the limit density functions (thick line) and the probability distribu- 
tions (thin line) as the initial state |yo(0)) = r [l,0] with = n/3. 



5. SUMMARY 

In the final section, we conclude and discuss the probability distributions of our walks. 
For the usual 2-state walk defined by U, the limit probability distribution of X x /x as % — > 
c° has a simple density function. On the other hand, if another matrix H operates the walk 
in certain intervals, the limit density function can be described by a 5-function and a 
combination of density functions. We found out a 2-state model whose limit distributions 
had the similar structures to the multi-state walks with more than four states defined in 
[17, 18, 19]. In Sect. 4, we presented the limit densities for some cases. Both coefficients 
Ai and A2 of the 5-function given by Eqs. (24) and (32) respectively were particularly 
independent from the initial state | Vb(0)) = T [ a iP]- One of the further problems for 
our model is calculation of the limit distribution for general m,n (m = 3, 4, . . .) case and 
m,n — > 00 case. 



APPENDIX 



A. CALCULATION OF THE LIMIT DISTRIBUTION 



In this appendix we calculate the limit distribution for m = 2,n = 3,4, . . . case. For 
the other cases, the distributions can be obtained similarly. Our approach is based on 
the Fourier analysis applied to QWs by Grimmett et al. [20]. For t = (n + 2)x + 2, we 
concentrate on the characteristic function E(e lzX '^ n+2 ^' c ) as z — > °°, where E(X) denotes 
the expected value of X. At first, the eigenvalues Xj{k) ( j G { 1 , 2}) of U (k) are computed 



as 



Xj(k) = -(-l)Vl - c 2 sin 2 k + ic sink. 



(37) 



The normalized eigenvector \vj(k)} corresponding to Xj{k) is 



\vj(k)) 



a/ 1 — c 2 sin 2 fc— (— l)Jccosk 
\ 2s 2 \ / l-c 2 sin 2 k 



se 



ik 



— ( — l) J \/l — c 2 sin 2 k — ccosk 



■ (38) 



By using the eigenvalues and eigenvectors, Eq. (20) is written as 



%(k))=(u(k)r(mu(ky) 2 \%(k)) 

={A 1 (fc)(' 1+2 ) T a 1 ,oAn-Ai(fc) ( "- 2)T a2,oAlTA 12 
+ (-l) T A 1 (fc)' lT (a 1 , A 12 + a 2 ,oA 11 )A 12 }|v 1 (fc)) 
+ {A 2 (fc)(' 3+2 ) T a 2 ,oATT 2 + A 2 (fc)(' 1 - 2 )^ 1 ,oA 11 A 12 
+ (-l) T A 2 (fc)" T (a 2 ,oA 12 -ai, AlT)A 12 } \v 2 (k)} , (39) 



where a jfi = a jfi (k) = (vj(k)\%(k)) mdA jlh =A hh {k) = (v h (k)\H(k)\v h (k)) (71,72 G 
{1,2}). Remark that \%(k)) = | ¥b(0)) = r [a,/3] (see Eq. (11)). The r-th moment of X t 
is expressed as follows: 



dk 
2k 



E(X[) = jyP(X f =x) = f (%(k)\ {D r \%(k))) 

= ((n + 2) r) r £ { Ifl^oA?! | 2 /*i (*)' + |a 2 ,oA^| 2 % (*) r } ^ 
+ {{n-2)x) r j ^{|«2,oAiiAi2| 2 ^iW"+|ai,oAiiA 12 | 2 /i2W)2^ 
+ («T) r y ^||(ai ; oAi2 + a 2 ,oAii)Ai2| 2 /zi(fc) r 

+ | («2,oAi2 - «i,oAiT) Ai 2 | 2 /i2 W"}^ +o(T r ), (40) 

withfy(Jfc) =DXj(k)/Xj(k),D = i(d/dk) and = f(f — 1) x • • • x (t - r+ 1). Noting 
that 



*y(*)=(-i)' 7 CC °f. 2 , (41) 

VI -c 2 shrfc 

4 //^ , /,x ccos 2 fc + z'srnfc\/l -c 2 sin 2 fc 

An(fc) =-A 22 (k) = , (42) 

VI -c 2 shr 

A12W =A 2 i(*) = , |5|C ° Sfc 9 , (43) 
Vl -c 2 shrfc 



in a similar way to Machida [21], we see that 



lim E 



(n + 2)x 



{ |«i, An I 2 ^i + l^.oA^ | 2 ^ 2 (/:) r ) 

71 ^ J 

+ J^\a 2 ,oATiA n \ 2 (^-^h^k) 
+ |ai ;0 AiiAi 2 | 



dk 
2n 



+ 



J ^||(ai,oAi2 + a 2 ,oAiiMi2| 2 ^^2*1 (*) 



+ I (02,0^12-01,0^11)^12!' 



iy{ 



n + 2 v ' J 2k 



n -\-2 n + 2 

f K (x;c)M 2 (x) + f K ( x; c)B 2 ( x; n) 

n n 

+ f K C-^x;c)B 3 ( r ^-x;n)}dx, (44) 
n — I n — I 



where 



/k(x;c) 



M 2 {x) 



k{ \ —x 2 )Vc 2 —x 2 



7 (-|c|,|c|)W> 



l-Ua\ 2 -\P\ 2 + W + a P ) s 



B 2 (x;n) = 



n + 2s 2 



n 



1+ l«| 2 -|/3| 



2 (a/3 + a/3>(3s 2 +l) 



• x 



2 |fl|2 , (a[5 + ap)s { s z 3 4(a/3 + a/3)s 3 x 



B 3 (x;n) 



-4j|a|--|/3| J f 
n + 2 



l-x 2 



n-2 



1 12 1012 (aj5 + aj5)s 
l + t |a| - \p\ l + ^ J - — 



1 - 

c z V c z 



(45) 
(46) 



(47) 
(48) 



and 7a (x) = 1 if x e A, 7a (x) = if x ^ A. By Eq. (44), we can compute the characteristic 
function £ , (y zX </("+ 2 ) T ) as t — > °°. Thus the calculation of the density function for 



m = 2, n = 3, 4, . . . case is completed. 
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